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Abstract 


An atmospheric attenuation model for the ultraviolet, visible, and infrared 
was developed in 1964, based on scattering (molecules and aerosols) and ozone ab- 
sorption. Since then more measurements have been made and our knowledge of 
aerosol attenuation has widened, These circumstances result in attenuation model 
changes which are relatively unimportant for most exploratory calculations. They 
can be significant, however, for long slant-path high<altitude applications entailing 
large zenith angles, factors which characterize, for example, the measurement 
geometries of rockets and satellites, Accordingly, a revision of the 1964 Attenua- 
tion Model is warranted. 

In this paper the optical parameters are computed spectrally and with altitude 
as follows: (1) pure air attenuation parameters are determined by utilizing Rayleigh 
scattering cross sections with molecular number densities from the standard at- 
mosphere; (2) ozone absorption parameters are derived based on Vigroux's coeffi- 
cients applied to a representative atmospheric ozone distribution; (3) seven sets 
of aerosol measurements are compared and a profile of aerosol attenuation coeffi- 
cients vs altitude is developed. Attenuation coefficients and optical thickness due 
to molecular, aerosol, and ozone attenuation are computed and tabulated individual- 
lv so that the influence of each can be compared. The newly derived tabulations 
permit various exploratory calculations, including horizontal, vertical, and slant- 
path transmission at kilometer intervals to an altitude of 50 km, individually for 
each attenuating component or for overall atinospheric extinction (molecular + 
ozone + aerosol). 
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Symbols ae 
A, Vigroux ozone absorption coefficient (em!) i 
Ds Ozone concentration (em /km) ‘ 
d Horizontal path length (km) . | 
H,, Aerosol scale height (km) [ 
h Altitude (km) m4 
K Mie scattering efficiency i 5 
m Aerosol index of refraction bg 
m, Index of refraction at sea level, air at 1013 mb and 15°C 2 . 
N, Ozone number density (cin) . a 
N, Aerosol number density (cm™°) ri 
N. Molecular number density (em™¥) . 
N, Molecular number density at sea level (em79) i. 
r Particle radius (microns) , 
th? Turbidity (6 ,/6,) iY 
T Temperature °K 
Th Horizontal transmission 
To-h Transmission between sea level and altitude h 
2 ae Transmission between h and space 
T Ah Transmission between two altitudes above sea level 
Ba Atmospheric ozone absorption coefficient (km™') 
By Aerosol attenuation coefficient (km7!) 
Bil ay) Mean of 79 profiles for A, = 0.55 microns (km"!) 
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4. Rayleigh attenuation coefficient (km!) 

Bost Extinction coefficient (km**) 

6 Depolarization factor ; 

6 Zenith angle : 

r Wavelength (microns or em) : 

on Wavelength 0.55 microns \ 

0. _ Rayleigh scattering cross section (cm?) : 

3 Ozone optical thickness from sea level to altitude h (O-h) 

3 Ozone optical thickness from altitude h tc space (h-« ) 

» Aerosol optical thickness from sea level to altitude h, (0-h) 
i Aerosol optical thickness from altitude h to <pace, (h- oo ) 

a Rayleigh optical thickness from sea level to altitude h, (Q-h) 
: Rayleigh optical thickness from altitude h to space, (h-© ) 


Sek Extinction optical thickness (molecular + ozone + aerosol) 
from sea level to altitude h, (O~h) 


Wie Extinction optical thickness (molecular + ozone + aerosol) 
from altitude h to infinity, (h-°:) 


y Aerosol size distribution function . 
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UV, Visible, and IR Attenuation for Altitudes to 
50 km, 1968" 
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1. INTRODUCTION 


In 1964, an atmospheric attenuation model was published (Elterman, 1964) 
which has been useful for a variety of exploratory calculations. Now a revision is 
warranted lor reasons given in the abstract and Section 4 of this report. In this 
revision most of the earlier material is presented in summary form. The section 
on attenuation by Rayleigh scattering, however, is retained because the content 
leading to the derivation of the Rayleigh parameters is useful. In one instance, 
due to existing interest, the material is expanded, i.e., the tabulations which com- 
prise the attenuation model now include aerosol and ozone optical thickness so that 
a comparison can be made of the relative importance of each attenuating component 
for vertical and slant paths. . 

The shortest wavelength considered is 0.27 microns. ‘The use of shorter wave- 
lengths would require a treatment of ©, absorption. Also, attenuation is suffi- 
ciently severe so that interest in the shorter wavelength region for purposes of 
ultraviolet transmission below 50 km probably is limited. The longest wavelength 
used is 4.0 microns. Calculations for longer wavelengths are complicated by the 
presence of absorption bands of 11,0, CO,, and their wings. In between, a total 
of 22 wavelengths is chosen (Table 1) within the atmospheric windows and for the 


ultraviolet region where ozone absorption is important. 


- (Received for publication 25 March 1968) 
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Conceptually, the attenuation model starts with molecular densities from the 
latest published U.S. Standard Atmosphere (1962) followed by the addition of ozone 
and ae€rosol components. The meteorological range (M.R.) at sea level corresponds 
to about 25 km at 0.554 wavelength. This choice serves a useful function because 
il permits including some important measurements conducted at\ = 0.55u. In 
addition, this wavelength customarily represents the phototopic region. 
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Table 1. Model Parameters as a Function of Wavelength 


cert RS tee ae le tla = Sain be moter Ae 


=26 eS 
1.00029668 | 8.960x 10 2.10 X 10 
1,00029475 | 7.64¢x1072° | 1.06x 102 
1.00029156 | 5.677x 10728 | 4.01 x 104 
1.00028902 | 4.310x1072° | 8.92x 1078 
1.00028609 | 3.334x1072° | 6,40x 1072 
1.00028531 | 2.622x10728 | 1,80x 1079 
100028392 | 2.091x10725 | o 
1.00028275 | t.e89x10725 | o 
1,00028052 | 1.038x 10725 | 3,50x 1072 
1.00027808 | 6.735x1072" | 3,48x 1072 
1.00027782 | 4.563x 10727 | 9,.20x 107? 
[ 1.00027607 | 3,202x1072” | 1.32x 1074 
i 1.00027630 | 2314x1072? | 6.20% 1072 
zi 100027878 | 1.714x10727 |. 2,30x 1072 
a 1,00027803 | 9,990x10°28 | 1,00x 107? 
100027451 | 6.213x 10728 
ce 1,00027387 | 3,216x 10722 
es 1.00027357 | 1.606x 10728 
bs 1.00027315 | 5.189x 10728 
: 1.00027282 | 1.817x 10728 
- 1,00027272 | 2,681x 10720 
1.00027269 | 1.571 x 1072? 


\ 7 Wavelength 
m, ~ Index of refraction (1013 mb, 15°C) 
o> Rayleigh scattering cross section 
A 
v 


~ Absorption coefficient after Vigroux for pure 03, smoothed 
values (1013 mb, 18°C) 
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2, RAYLEIGH PARAMETERS 


A fundamental requirement for generating an accurate series of Rayleigh 
parameters is an exact determination of the index of refraction for the wavelengths 
of interest, With this known, the Rayleigh cross sections can be computed. This 
in turn permits computation of the Rayleigh attenuation coefficient and its vartation 
with altitude as well as corresponding optical thickness values, 

The index of refraction for a standard atmosphere (1013 mb, 15°C) specifically 

for the desired wavelengths used is determined by Edlen's (1953) expression 


(m, = 1)107% = 6492.8 + ~2240810, , 25540, (1) 
146.6 (0%) 4a - (074) 


where m, ? refractive index , 
4 = wavelength (microns). 


Penndorf's (1957) computations using Eq. (1) demonstrate that the effect of water 
vapor oan be neglected andthe derived m, values have negligible error for the 
spectral range from 0.2 to 20.0 microns, 

The Rayleigh cross section is expressed by 


3.2 4.2 
i (a) a 82 (m5 1) : 6 + 36 (2) 
r 3* n2 6 = 7 


o, = the Rayleigh scattering cross section (em?), 
A = the wavelength (cm), 
m, »® the index of refraction of air at 15°C and 1013 mb pressure, 


» the molecular number density at sea level for a standard 
atmosphere (em=3), = ° 


6 = the depolarization factor, 


The term (6 + 38)/(6 ~ 78) accounts for the degree of depolarization attributable to 
the anisotropy of the atmospheric molecule. The depolarization factor has heen 
determined by-calculation and by laboratory measurement. ‘I'he latest work of 
Gueker and Basu (1953) ylelds 620.035. The wavelengths of intereal with the 
indices of refraction and Rayleigh cross sections [computed from Iqs, (1) and (2)) 
are listed in Table 1, and plotted in Figures t and 2. 

Using the scattering cross sections, the Rayleigh coefficients are obtained 
with 
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Figure 1. Index of Refraction for 1013 mb and 15°C (Table 1), 
cooo Represent Attenuation Model Wavelengths 
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“ Figure 2, Rayleigh Cross Section vs Wavelength (Table 1), 
ay 0000 Represent Attenuation Model Wavelengths 
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+ BRAY = OLA) © NOChp (10° cm/km) , (3) 


where 


Bo. = the Rayleigh attenuation coefficient (km” F 
a. the Rayleigh scattering.cross section (cm?), 
N, * the molecular number density (om73) . 


The values of N Ch) needed for Eq. (3) were obtained from the U.S. Standard 
Atmosphere and are listed in Table 2. This expression is used to compute an array 
of Rayleigh attenuation coefficients as a function of altitude for each wavelength. 

With the Rayleigh attenuation coefficients determined, the optical thicknesses 
from sea level to some altitude h are computed by 


h 
ThA) = Di kh) Oh ; (4) 
0 


t, = Rayleigh optical thickness (0~h), 


B. = mean of the Rayleigh attenuation coefficients (km~') 
for each altitude increment, 


4h = altitude increment chosen as one km for these calculations. 


The Rayleigh optical thickness for altitude h to space is obtained by the rela- 
tionship , : 

mich, A) = te »A) = 7 (ha) sete: : i (5) 
where : ; a: % 
71.(h) s Rayleigh optical thickness (h-e) , 
t AC «) = Rayleigh optical thickness (0-o)' . - eS \ 


The term Flee) was obtained by using Eq. (4) with the limits set between 0\and 
80 km. Above 80 km, Stergis' (1968) calculations, based onN,,O,, and O \as 
the principal atmospheric constituents, yield 


= 3.6X107° dX 50.4 

J 5.gnayen : 6. TX 1077 = 0.6 ; 

80 21X10 ° ,dA 20.8p 
These values approximate a constant 1075, that of the Rayleigh optical thickness 
for unity air mass. For our purposes then the integral can be neglected because 
the constant is small and applies to all waveleng'"s of interest. 


Table 2. Model Parameters ag a Function of Altitude 


2.547 x 10"? 3.56 X 107 


2. 311 4.26 gins 2 _ , i 
2:093 2.03 ° ¢ i 
1.891 2.50 ! aur 
1.704 2.26 e 3 
1,531 2.210 “ore 3 i 
1.373 2.16 . : 
1.227 2.23 
1.092 4, 2.28 ‘ 
9.712 10 2.81 
8,598 3.50 
7,585 4,60 
6,486 r 6.21 
5,543 8.45 ; , 
4.738 9.57 DAN { 3 
4,049 9.94 “2 aur : 
: 3.461 1.03 X 10 ; 
2.959 1.1 i 
2.162 1.42 
é 1.849 1.64 
1.574 1.84 | 
1.341 1.97. 
1.144 49 1,98: { 
9.760 X 10 1,93 
8,335 1.80 | 
e 7.123 1.63 
fe 6.092 1.41 } 
he 5.214 1.23 ' 
a 4,466 1.07 4 
a 3,828 9.03 X10 
| 3,283 7.93 
ao 2.818 6.82 i 
iF 2.406 5.82 H 
it. 2.066 4,85 } 
1.760 4,31 
1.509 3.61 | 
H 1,296 3.02 i 
1.118, 1¢ 2.83 
i 9.620 X 10 2.17 | 
f * 8, 308 1.86 
7.187 1.52 . 
6,227 1.19 “4 Vat oe . 
5.404 9.30 X 10 
4,697 7.44 
4,088 5,76 ; 
3,584 4,46 } 
3.112 3,53 : 
2.738 2.79 
2.418 2.23 
2.135 1.86 


h= Altitude; Np— Molecular number density; 
D3- Ozone equivalent thickness 
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3, ABSORPTION PARAMETERS FOR ATMOSPHERIC OZONE 


This section is a summary of the material used in the 1964 Attenuation Model, 
The parameter for determining Og apsorption ag a function of altitude ia -the 
atmospheric ozone absorption coefficient expressed by: , 


Beha) = ALA) Dah), (8) 


B 3 = atmospheric ozone absorption coefficient (km7!), 
A. = the pure ozone absorption coefficient (em"}) after Vigroux, 


Vv 
D the ozone equivalent thickness (cm/km). 


3 s 
Thus, the Vigroux coefficients (1953) listed in Table 1 in conjunction with the ozone 
concentrations, Table 2 and Figure 3, permit the computation of an array of atmos- 
pheric ozone absorption coefficients to 50 km for each of the desired wavelengths. MS 
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Figure 3, Representative Atmospheric Ozone 

Concentration Profile (Table 2). Values for: 

0 to 30 km, based on Handvook of Geophysics 
and Space Environment (1965) 

30 to 40 km, interpolated 

40 to 50 km, based on London, Ooyama, and 
Prabhakara (1962) : 
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The ozone optical thickness from sea level to altitude h, t gl h,A), and the 
ozone optical thickness from some altitude h to space, rig(hid), are included in : 
the model tabulations, ‘‘he expressions for deriving these parameters have the 
same form as Eqs. (4) and (5). 


4, AEROSOL ATTENUATION 


Of the various methods used to investigate aerosol attenuation, for the present 
we will consider optical techniques only because they are suited to this type of study, 
In this country, Newkirk and Eddy (1964) used solar aureole photometry; Penndorf 
(1954) analyzed solar radiation measurements from aircraft altitudes; Elterman's 
results (1966a,b), with searchlight probing, comprise a substantial number of pro- 
files acquired in New Mexico for altitudes to about 34 km, In Australia, Crosby 
and Koerber (1962) used a balloon=borne integrating nephelometer. In the U.S.S.R., 
Kondratiev, et al. (1987), conducted balioon solar transmission measurements ; 
Feoktistov (1965) analyzed photographs of the earth's horizon from the spacecraft 
Voschod; Rozenberg, et al. (1980) (1966), acquired their results with searchlight 
probing. The various results, as shown in Figure 4, were made comparable af. 

“4 = 0, 55h by using the empirical relationship that the aerosol attenuation codffi-~ 
cient is inversely. propbrtidnal to wevesrae For reasons of clarity, a substantial 
body of resulta wat ‘Not included in Figure 4, ‘as for example the twikight ‘measure- 
ments by Rozenberg (1965), Volz and Goody (1962), the searchlight measurements 
by Spankuch (1967), analysis of twilight abreole photographs from. the spacecraft 
Vostok-6 by Driving (1966), the aircraft measurements of sky brightness by San- 
domirski, Al'tovakaia and Trifonova (1964), and aircraft nephelometry by Waldram 
(1945), Also, interesting results in the form of relative values have been obtained 
with optical techniques: the twilight measurements by Bigg (1964), the laser beam 


_ backscatter by Collis and Ligda (1966), by Paeeree Kent, and Wright (1967), 


and by Grams and Fiocco (1967), 

A consideration of all results shows, \as does Figure 4, that the aerosol atten- 
uation coefficient is a strongly fluctuating parameter and that average values based 
on an adequate number of measurements are necessary in order to establish a repre- 
sentative profile. The recent searchlight probing measurements (Flterman, 1966a 
and 1966b) appear representative based on several considerations. First, each 
profile was acquired by continuous measurement through both troposphere and 
stratosphere and with an altitude resolution approximating one km. In addition, a 
total of 119 profiles comprising absolute values of aerosol attenuation coefficients 
were determined for various times throughout the year. This represents a sub- 
stantially larger sample than previously published, Further, such a quantity of 
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AEROSOL ATTENUATION COEFFICIENT (km~') 


107! 


10°? 


io’ 


, 35 
ALTITUDE (km) : 


e 4, Aerosol Attenuation Coefficients ve Altitude ath; = 0.55y. 

Conjparison of results: 

XxX/AX solar aureole, 2 balloon flights, Newkirk and Iddy (1964); 

44/44 solar radiation measured from aircraft, mean of 8 flights, 
Penndorf (1954); 

searchlight probing, mean of 105 profiles, Flterman 
(1966a) and (1966b); 

balloon integrating yma mean of 14 flights, Crosby 
and Koerber (1962); 

solar radiation measured from balloons, mean of 3 flights, 
Kondratiev et al. (1967); 

spacecraft horizon photography, analysis Sf 4 frames, 
Rozenberg (1966), Feoktistov (1965); 

searchlight probing, mean for 5 nights, Rozenberg (1966) 
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results readily permits statistical treatment. KFinally, we note that the mean of the 
119 profiles falls reasonably well within the values determined by other researchers, 
a circumstance which tends to provide a measure of comfort. 

In considering the suitability of the 119 profiles, -extensive averaging is re- cae 
quired and this tends to wash out features easily noted in the individual profile. 
We present, therefore, in Figure 5 a single profile, chosen because its properties 
are readily evident and also because it is similar to the overall average. The fea- 
tures can be made to emerge more prominently if the aerosol coefficienta are used 
to compute a turbidity profile, t (hry) +B (hd )/6,(h, a4) , where By and B. are 
the aerosol and Rayleigh coefficients respectively and ST = 0,55u (Figure 6). 

Since volcanic dust in the atmosphere can have a residence time of several 
yeas, the effects of the Mt. Agung eruption (March 1963) must be considered. 
‘The direct measurements of Junge, Chagnon, and Manson (1961), Friend (1965), 
Mossop (1964), and Rosen (1968), collectively considered, before and after this 
event, show evidence of change in the stratospheric aerosol content. The observa- 
tions of the twilight sky by Volz (1965) and Meine) and Meinel (1964) also show 
augmentation of stratospheric particulates. Since the searchlight probing measure- 
ments yielded absolute values of aerosol attenuation coefficients, the most suitable 
parameter to use for examining this feature quantitatively is tne stratospheric Tie i . 
aerosol optical thickness for the altitude region between the tropopause and 25 km, . 
The reason. for choosing the latter altitude limit will be discussed later, Accord- 
ingly, all profiles were placed in time-sequential groups determined by the simi- — 
larity of the stratospheric ‘dust feature. Then the optical thickness was computed by 


h 


n 
=i YY suman , (7) 
ial hy, 


where n is the number of profiles in the group, hy is the altitude of the tropopause, 
ho the 25 km altitude, §' the mean aerosol attenuation coefficient (within each pro- 
file) for the altitude interval, and Ah the altitude intervals used for computing 

the profiles. The results of this computation are presented in Table 3, The tabu- 
lation demonstrates a relatively high level of stratospheric dust for the December 
1963 to March 1964 period. Beginning with April 1964, dust abatement and a gen- 
erally stabilized level are in evidence. The mean optical thickness of Group (B+C+D) 
is 26 percent less than that of Group A. Since Group A entails a period of abnorm- 
ally high aerosol content, its profiles are not representative. These results are 

in satisfactory agreement with the findings from the direct measurements of the 
authors mentioned. 
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Figure 5, Single Profile Bph, At) for 11 April 1964 at 02:00 MST, 
Similar to Mean of 79 Profiles, \1 = 0.55y. 
Surface to 5 km — convective region; 
\ 5 to 11.7 km — troposphere dust layer; 
11.7 to 23,8 km-— stratosphere dust layer; 
25.6 km — upper altitude maximum : 
+++ + aerosol (measurements) 
~ molecular (computed) 


i 
[ 


ALTITUDE (km) 


12 


40 


25 


2 3 
TURBIDITY Bp/B, 


Figure 6, Single Turbidity Profile 8 (h, A978 (hay) 


for lt April at 0200 MST, 
+++++ Measurements 


‘(See caption pertaining to Figure 5) 
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Table 3, Aerosol Optical Thickness as a Measure of Volcanic Dust, d= 0. 554 


Period Number of Mean Optical Thickness 
(Inclusive) Profiles (approx, 12-25 km) 


1963 — Mar 1964 
1964— Sep 1064 
1964 = Nov 1964 
1964 — Apr 1965 
1964 — Apr 1965 


The considerations discussed above justify the selection of the Group cose 
profiles as a basis for developing representative aerosol attenuation parameters. 
It will be convenient to designate the Tz profile average for AY 20.56 as 
B (rs Ay) (Figure 7). This profile can be extended to encompass a larger altitude 
range by using the scale height relationship, 


As - * (hy = hy) 
By(hos Ay) * B (hy Ay) exp | - ree? ‘ (8) 


Penndorf'a study (1954) shows that for the loweat 5 km, the aerosol coefficients fall 
off exponentially with a scale height 0.97<H <1.4 km. We resort to the use of 
his mean value, Hp ‘21.2 km, to extend the 6 (th, a) profile from 3.7 km to sea 
level. Thig results in aerosol coefficients which are identical to those of the 1964 
Attenuation Model for altitudes 0 to 3 km (Table 4,11). 

Above the convective region, up to the tropopause the aerosol coefficients show 
a moderate gradient which is in close agreement with Penndorf's analysis (1954) of 


- the vertical distribution of aerosols in the troposphere. Also, this section of the 


profile is based on high signal to noise measurements and extensive averaging, 
factors which contribute to its reliability, Additionally, this altitude region, being 
above the convective layer, is ciaracterized by aerogol conditions which tend to be 
independent of surface terrain, These considerations suggest that the shape and 
values of the 8 (h, 1) profile for this altitude region are realistic. 

Above the tropopause, up to 25 km, the cocfficients are larger than those de- 
rived for the 1964 Attenuation Model by a factor of about 20 (at 20 km for AY = 0, 55m). 
This difference may be attributed in part to the intrinsic difficulties of converting a 
size distribution to an optical parameter, that is, establishing the radii mits, - 
particle shape, chemistry, and index of refraction. Relitive to profile shape, a 
turbidity maximum dominates at about 19 km (Migure 8). ‘The measurements of 
Rosen (1968) and those of Volz (1968) for 1964 to 1968 are in satisfactory agreement 
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Figure §, Mean Turbidity Profile, 4 (hy) 7/8 (hed 4) 
(See caption for Figure 7) P 
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with the mean optical thickness and shape of the Group (B+C+D) profiles, Table 3. 
Also their results show that the stratospheric dust level has remained approximate- 
ly constant from 1964 to the time of this writing. The B (hay) curve (Figure 7) 
shows this dust feature as the knee of the profile rather than a massive layer fea~- 
ture indicaied by the turbidity profile. This conceptual relationship suggests that 
over-emphasis is pussible when dealing only with turbidity profiles. aaa 

The turbidity profile shows an upper stratospheric maximum with its lower 
terminus at 25 km. ‘This altitude then, was the basis for choosing the upper strato- 
spheric dust limit in Mq. (7), The maximum at 26 km occurs with sufricient fre- ; 
quency to be easily identified in the turbidity profile. The existence of such an i 
aerosol concentration above 20 km is supported by the analysis of satellite photog- 
raphy reported by-Mateer, Dave, Dunkelman, and Evans (1967). 

To establish upper altitude aerosol coefficients, a least square fit was com- 
puted for Bplhady) from 26 to 32 km (Figure 9). The result, Hi, = 3.75 km (in 
effect derived from 790 measurement points), was used in Eq, (8) to extend the 
values to 50 km. Miller (1967) obtained Hy = 3.25 km from a thorough analysis 
of rocket measurements acquired in 1964 for this altitude region. The overall 
profile from sea level to 50 km is presented in Figure 10 (Table 4.11). Since the 
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Figure 9, Expanded Scale for 26 to 32 km Altitude Region 
Showing Hp = 3.75 km for 79 Profile Mean; Aerosol 
Attenuation Coefficients, Bplh,A4) va Altitude; Least 
Square [it Used to Extrapolate to 50 km; Ay = 0.554u 
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Figure 10. Cornparison of Profiles. Aerosol attenuation coefficients 
6,(h, 4), with molecular 6y(h, Ay). The 1964 profile shows interpola- 
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turbidity is proportional to the mixing ratio, the diminishing values for the extrapo- 
juiiun imply the aciveul puuve caisia below OO km. Should the source be of meteoric 
origin, tie mixing ratio would tend to be constant or increase for altitudes 30 to 50 
km. 

Thus far we have selected a set of measurements for Le 0.55y and provided 
reasons for ‘ts use. It would be in order now to examine some expressions leading 
to corresponding aerosol profiles for other wavelengths. If we consider a real i 
atmosphere, the aerosol sizes within unit volume can be described by a size dis- a 
tribution function, ¥(r). Various size distribution functions are in use: the Junge ‘ 
type power law (1963) with a choice of exponents discussed in detail by Bullrich 
(1964), a similar distribution modified by gaps observed by Fenn (1964), a log- 
Gaussian distribution used by Foitzik (1965), a composite distribution with com- 
ponents from several types. The optical-particle size relationship utilizes (r) 


such that 
Ne : 
sym) = fi otmra) aNgind , () | 
Neg i 
aN, = Ny wr)dr , (10) ; 
Noh) = ONAH) as 


B.. is the aerosol attenuation coefficient, m is the index of refraction, N,, and 
Nrp are the aerosol number density limits established by the radii limits r, and 
rg, 9, is the aerosol cross section for each particle, N, is the total number of 
particles between r; and rp. For-a given altitude, No actually is proportional 
to the particle number density between ry and Pos Binge the same size distribu- 
tion function applies to all altitudes (an assumption), Eqs. (9), (10), and (11) are 
combined 


ro 


B(h, d) = CNA(h) Sf on(red) wir) dr. (12) 


ry 


Here, CN.(h) is placed outside the integral which now contains only factors that 
are independent of altitude; also, m is removed because subsequent considera- 

‘ tions will pertain to particles without any distinction in refractive index. If 
Eq. (12) is normalized to sea level conditions, the integral cancels out. Then 
generally for the various wavelengths A, and specifically for Ay =0,.55p , we have 
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B (h,a) B_(h. 2.) N (hd 


BOAT * BOX) * RO (13) 


or 
6(0,%) 
Athy») = BOR Ath, A). (14) 


Equation (13) has been derived in this manner to demonStrate its compatibility with 
particle size considerations. Sea level conditions have been researched extensively 
by Curcio and Durbin (1959), Curcio, Knestrick, and Cosden (1961), Knestrick, 
Cosden, and Curcio (1961), Dunkelman (1952), Baum and Dunkelman (1955), The 
6(0,4) values fora 25 km M.R., based on the results of these authors,are shown 
in Figure 11 (see also Elterman, 1964). Utilizing these results, in conjunction 
with the Ath, AD profile (Table 4.11), all requirements for the right-hand side 

of Eq. (14) are satisfied and an array of aerosol attenuation coefficients can be 
computed for all altitudes and wavelengths of interest. 

The aerosol optical thickness from sea level to altitude h, r.(h,2), and the 
aerosol optical thickness from some altitude h to space, 7'(h, Ni » are included in 
the model tabulations. The expressions for deriving these parameters have the 
same form as Eqs. (4) and (5), 
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Figure 11. Aerosol Attenuation Coefficients §,(0,A) vs 
Wavelength at Sea Level for a Meteorological Range 
Approximating 25 km. 

A -— derived from Baum and Dunkelman (1955) . 

B ~ contained in Curcio, Knestrick, and Cosden (1961) 


20 
3. ATMOSPHERIC EXTINCTION 


In this section, three sets of extinction parameters are considered: extinction 
coefficient, extinction optical thickness from sea level to a desired altitude, and 
extinction optical thickness from a desired altitude to space. 

The atmospheric extinction coefficient 6 ext is the sum of all the attenuating 
components: 


Bexgt(Bed) = B(h,A) + fglh, A) + Bo(h,Ad (18) 


The extinction optical thickness from sea level to altitude h, Text itr) » and 
the extinction. optical thickness from some altitude h to space, Th (hn) , are in- 
cluded in the tabulations of the attenuation model, The expression for deriving 


these parameters has the same form as Eqs. (4) and (5). 


6. EXPLORATORY TRANSMISSION CALCULATIONS 


Using the derived tabulations that follow, some exploratory calculations with 
extinction parameters (for any of the wavelengths) are,démonstrated. Rayleigh, 
aerosol, and ozone parameters can be used similarly. 

For horizontal transmission (T,) over a path (d) at any altitude (h), the ex- 
tinction coefficient 


T, = exp [ = Bexatnadea | : (16) 


For vertical and slant=path transmission from sea level to a given altitude, at 
zenith angle @ for all wavelengths of interest 


Tyan = OxP [- Toxtih) * Bec | . (17) 
For vertical and slant=-path transmission between two altitudes above sea level 
T an = @XP - [ rent(ltg) - rext(hy) | *sec@ , (18) 


For vertical and slant~path transmission from a given altitude out into space 


Theo 7 XP [- 7 hy (h) sec a] . (19) 
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7% CONCLUDING REMARKS 


The procedure for developing the aerosol attenuation profile is summarized as 
follows: a 

(1) Various studies were compared and of these a set of measurements 
selected. 

(2) The choice of measurements (comprising 119 profiles from 2.76 to 34,4 km) : 
was examined statistically. This resulted in the elimination of 40 profiles 
(December 1963 to March 1964 inclusive) characterized by a high volcanic 
dust component. 

(3) The mean of the 79 remaining profiles ‘was extended to sea level and to 
50 km respectively by reasonably supported extrapolations. : 

(4) The overall profile then was developed laterally to obtain 21 additional 
profiles for the wavelengths of interest. 

A significant aspect of the procedure is that the wavelength-height array of param- 
-eters was derived independently of the assumptions associated with conversion of 
a size distribution to an optical parameter. . 

For most purposes, calculations using the new parameters will be affected 
only moderately. For example at Ay 2 0,55p, the 1964 Attenuation Model provides ; 
an extinction optical thickness, tee 0.331, for a vertical air mass. The new ae 
parameters yield Tiegh = 0.379, resulting in a transmission change of about 3-1/2 
percent. However, for long path horizontal transmission calculations above 5 km 
and for long slant-path calculations entailing large zenith angles, the new aerosol 
parameters can function more significantly. 

As mentioned previously, the Rayleigh and ozone parameters are unchanged. 
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8 TABULATION OF PARAMETERS 


Tables 4,1 to 4,22, which follow, comprise the atmospheric attenuation model, 
Exponents are in computer notation; for example, read 2,86-3 = 2,86 X 1073 and 
2.86 3.8 2.86 10°, 
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